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ABSTRACT 

We report on the analysis of high-resolution infrared spectra of the newly discovered 
brown dwarf e Indi B. This is the closest known brown dwarf to the solar system, 
with a distance of 3.626 pc. Spectra covering the ranges of A2. 308-2. 317/xm and 
Al. 553-1. 559/im were observed at a resolution of A/AA=R=50,000. The physical 
parameters of effective temperature and surface gravity are derived for e Ind Ba by 
comparison with model spectra calculated from atmospheres computed using unified 
cloudy models. The results are Teff= 1500±100K, log g= 5.2ib0.3 (in units of cm 
s~^), placing it in the critical boundary between the late-L and early-T dwarfs. The 
high spectral resolution also allows us to measure an accurate projected rotational 
velocity, with vsin(z.)= 28±3 km s"-*^. Combined with a published luminosity for e Ind 
Ba (with log(L/L0)= -4.67), the derived parameters result in a "spectroscopic" mass 
estimate of ~32Mjupitcr, a radius of ~ O.O7R0, and a maximum rotational period of 
~ 3.0 hours. A compilation and comparison of effective temperatures derived from 
spectroscopy using model atmospheres versus those derived from luminosities and 
theoretical Mboi-radius relations reveals a systematic disagreement in the TgS scale. 
The source of this disagreement is unknown. 

Subject headings: infrared: stars — stars: brown dwarfs, fundamental parameters, 
individual (e Ind B) 

1. INTRODUCTION 

The recently discovered nearby T-dwarf companion to e Indi (Scholz et al. 2003) will be 
important in improving our understanding of the behavior of the brown dwarfs that fall within 
the newly defined L and T spectral types (e.g., Kirkpatrick et al. 1999; Geballe et al. 2002). 
The T-dwarf reported by Scholz et al. (2003) was discovered by Volk et al. (2003) to be a close 
optical double consisting of an early T dwarf (e Indi B) and a late T dwarf (e Indi C) separated 
by 0.6 arcseconds. With an accurately known distance of 3.626ib0.009 pc, e Indi Ba and Bb are 
the nearest known brown dwarfs. They share a common proper motion with the K5V star e Indi, 
lying at a projected distance of ~1460 AU from their presumed primary star. The age of e Indi 
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itself has been estimated by Lachaume et al. (1999) to be ~0.8 to 2.0 Gyr: this age estimate is 
based upon its rotational velocity and Ca II K-line emission, e Indi Ba and Bb are thus brown 
dwarfs with very well-defined luminosities and approximate ages. 

We present the first high-resolution infrared spectroscopic observations of e Indi B. Synthetic 
spectra computed from unified cloudy models by Tsuji (2002) are compared to the observed 
high-resolution spectra; these comparisons are used to derive the stellar parameters of effective 
temperature (Tegf), surface gravity (defined as log g), and projected rotational velocity (vsin(i)). 

2. OBSERVATIONS 

High-resolution infrared (IR) spectra were obtained on e Indi Ba using the 8.1 m Gemini South 
reflecting telescope and the NO AO Phoenix spectrometer (Hinkle et al. 1998). This instrument is 
a cryogenically cooled echelle spectrograph that uses order separating filters to isolate individual 
echelle orders. The detector is a 1024x1024 InSb Aladdin II array. The size of the detector in the 
dispersion direction limits the wavelength coverage in a single exposure to about 0.5 % (1550 km 
s~^ or ~0.0120//m at 2.3/im and ~0.0080/xm at 1.6//m). One edge of the detector is blemished 
so the wavelength coverage on all but the brightest source is typically trimmed a few percent to 
avoid this area. The spectra discussed here were observed with the widest (0.35 arcsecond) slit 
resulting in a spectral resolution of R=A/AA= 50,000. Two spectral regions were observed, with 
one centered at A=2.314;tim and the other centered at 1.555//m. These spectral regions sample 
crucial diagnostic lines from the molecules CO, H2O, and CH4. 

At the time of our observations we were unaware of the existence of e Indi Bb. Following 
the discovery (Volk et al. 2003) we re-examined our acquisition images. Acquisition images of e 
Indi B at 1.647 /im taken on 2003 August 13 under good conditions (0.4" FWHM delivered image 
quality (DIQ)) confirm the Volk et al. (2003) detection. The companion can also be seen in 1.558 
)um images on 2002 December 29, although it is less well-resolved due to inferior seeing (0.8" 
arcsec DIQ). It is, however, barely perceptible at 2.321 /xm (2003 January 16) with 0.4" DIQ. The 
issue relevant to the current investigation is the extent to which the spectrum of e Indi B might 
suffer contamination from a nearby companion. The observed image profile at 1.647 /xm is well fit 
by a model in which the companion is fainter by 1.9 mag and is 0.65" away at a position angle 
of 125 degrees. At 2.321 fim. the magnitude difference is > 3 between the two stars, so that the 
companion could contribute no more than 6% to the spectrum in this region. Our data prevent 
an estimate at 1.558 jum but Volk et al. (2003) report a difference of 1.3 magnitudes. In addition 
the position angle is at 35 degrees to the 0.35 arcsecond slit. Detailed examination of the spectral 
images showed no trace of e Indi Bb at the location expected so we are confident that our spectra 
are contributed almost entirely by e Indi Ba. 

Each program star was observed along the slit at two or three separate positions separated 
by 4" to 5" on the sky: the delivered image FWHM at the spectrograph varied from 0.25"-0.80" 
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during the nights that spectra were taken, so stellar images at different positions on the slit were 
well separated on the detector. Equal integration times were used for a particular program star 
during a particular set of observations. With this observing strategy, sky and dark backgrounds 
are removed by subtracting one integration from another (the star being at different positions on 
the detector array). During each night, 10 flat-field and 10 dark images were recorded for each 
given wavelength setting of the echelle. A hot star, with no intrinsic spectral lines in the regions 
observed, was also observed each night in each observed wavelength region. 

Two-dimensional images were reduced to one-dimensional spectra using an optimal extraction 
algorithm described in Johns-Krull, Valenti, & Koresco (1999). Wavelengths and telluric 
corrections in the 2.31 jum region were determined by fitting observed telluric features with a 
scaled atmospheric transmission function from Wallace & Hinkle (2001). There are no significant 
telluric features in the 1.56 /nm region, so the wavelength scale was determined by matching a 
spectrum of HR 1629 (K4 III) with an IR FTS atlas of Arcturus (Hinkle, Wallace, & Livingston 
1995) 

Figure 1 illustrates the combined and reduced spectra for e Indi Ba in both the 2.313/um 
and 1.556/xm regions, with the wavelengths plotted as air wavelengths. The data points have 
been smoothed over the slit width of 4-pixels and the final signal-to-noise ratio is about 30-40. 
The 2.31/im region contains strong vibration-rotation lines from the first overtone bands of CO 
(here (2-0) lines from ^^C^'^^O), as well as some weak, blended H2O features, and weak absorption 
from methane. Detectable spectral features at 1.55/xm are not as strong as at 2.31/xm, nor as 
well-defined as the individual vibration-rotation CO lines, and consist mostly of blended H2O 
features. Two of the stronger features can be assigned to mainly two H2O lines (as marked in 
Figure 1) on the basis of the study of H2O by Tereszchuk et al. (2002). 

3. ANALYSIS AND DISCUSSION 

The observed spectra of e Ind Ba at 2.313/xm and 1.556//m are compared to synthetic spectra 
calculated from model atmospheres as discussed by Tsuji (2002). These models are so-called 
"unified cloudy models" in which dust is allowed to exist in the photosphere over a limited range 
defined by a condensation temperature, Tcond) and a critical temperature, Tcr, such that dust 
is found in the region of Tcr < T < Tcond- At the critical temperature, dust grains become so 
large that they precipitate from the the photosphere. The models employed in this analysis are 
computed with plane parallel geometry, in hydrostatic equilibrium, and have solar abundances. 

Figure 2 illustrates a comparison of observed and synthetic spectra, for both the 2.313)um 
(top panel) and 1.556jLtm (bottom panel) regions. The 2.313jLtm spectrum contains strong 
12q16q (2-0) lines and some weak H2O features. The 1.556;um spectrum exhibits primarily H2O 
absorption, with these features composed of many blended individual spectral lines. In the top 
panel, the comparison synthetic spectra span effective temperatures from 1400K to 1800K and 
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these models have surface gravities of log g= 5.5 (in units of cm s~^). This particular spectral 
region is illustrated as any CH4 absorption beginning near A~ 2.3158/xm, and clearly apparent in 
the models with Tcff=1400K or 1500K, is very temperature sensitive. Its observed absence (or 
extreme weakness) in e Ind Ba indicates that Tefif= 1600K; higher effective temperatures begin to 
produce CO lines that are too weak. The CO lines that dominate this region are not very sensitive 
to gravity over the expected range for the brown dwarfs with these approximate temperatures; 
however, the 1.556;um region features dominated by H2O are more sensitive to surface gravity, as 
shown in the bottom panel of Figure 2. Here, a slightly lower effective temperature is derived, with 
Tcfj=1400K, and the gravity sensitive H2O features indicate that log g~5.0 to 5.5. Higher effective 
temperatures produce H2O absorption features that are too weak for any reasonable surface 
gravity, while temperatures much lower than Teff=1400K (say 1300K) produce increasingly strong 
CH4 absorption, which looks nothing like the observed spectrum of e Ind B. This effect is shown 
by the Tcfj =1300K model spectrum in the bottom panel of Figure 2 that is offset vertically from 
the observed spectra and 1400K models. The offset is done as the different absorption features in 
the 1300K model (that are caused by increasing CH4 and decreasing H2O absorption), if overlaid 
on the observed spectrum, would produce merely confusion. 

Fits to the line profiles (primarily the strong CO lines), as shown in Figure 2, also yield 
the projected rotational velocity, with vsin(i)=28 km s~^ (with an uncertainty of ±3 km s~^). 
Wavelength shifts between observed and synthetic spectra also provide an accurate radial velocity 
for e Ind B, which we find to be Vheiiocentric= -41.0±0.7 km s~^. This radial velocity is very close 
to the published value for e Ind A's velocity of -39.6±0.8 km s^^ from Wielen et al. (1999). The 
similarity of radial velocities for both e Ind A and Ba strengthens their physical association, as 
argued by Scholz et al. (2003) based on their respective distances and proper motions. 

The combination of both the 2.314//m and 1.555/xm high-resolution spectra and their 
comparison to synthetic spectra result in values for temperature and gravity in e Ind Ba to be 
Teff =1400-1600K and log g=5.0-5.5. Taking the average of these values as being the best estimates 
we find an effective temperature of 1500K and a gravity of log g=5.25. Scholz et al. (2003) 
used photometry to derive the luminosity of e Ind Ba and found it to be log(L/L0)=-4.67. This 
luminosity can be combined with our estimates of temperature and gravity to yield a spectroscopic 
mass estimate of 32Mjupitcr- Scholz et al. (2003) derive a mass of 40-60Mjupitcr , which is close to 
the mass derived from the high-resolution spectra. In addition, given the luminosity and effective 
temperature, the radius of e Ind Ba can be estimated and then combined with the projected 
rotational velocity to yield a maximum rotational period. With log(L/L0)=-4.67 and Teff= 
1500K, a radius of (R/Ro)= 0.069 is derived. Given vsin(6)=28 km s~^, then the maximum 
rotational period for e Ind Ba will be 3.0 hours. 

The effective temperature derived by Scholz et al. (2003) is 1260K and results from the 
luminosity combined with the distance and a radius defined by a theoretical Mboi-radius relation 
that results from structural models of brown dwarfs. This Tcfj is somewhat lower than our value of 
1500K, however this difference is typical of the differences found to date between spectroscopically 
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derived effective temperatures and those derived from structural models. Tfie differences in Teff 
also indicate differences in the implied radii from the two methods. In our case, the higher 
spectroscopic Tcfj for e Ind Ba requires a smaller radius to support its luminosity. The two 
techniques are complementary, in the sense that the radius does not enter into the computation 
of the plane parallel atmosphere, but does in the Mboi-radius relation. Age uncertainties can 
also affect derived physical paprameters, especially for the structural models. Larger sets of 
comparison measurements of physical parameters derived from both model atmospheres and 
structural models will improve our understanding of the root cause of these differences. A 
comparison of the temperature scales is illustrated in Figure 3, where spectroscopic Tcfj's are 
plotted versus structural Tegf's; the spectroscopic temperatures are taken from Basri et al. (2000), 
Leggett et al. (2001), and Schweitzer et al. (2002), while the structural temperatures are those 
from Dahn et al. (2002). The Dahn et al. Teff's also use a Mboi-radius relation from structural 
models in deriving temperatures. Despite using different sets of model atmospheres (with differing 
treatments of dust), there is a clear trend in the differences between the spectroscopically derived 
effective temperatures when compared to the structural temperatures from Dahn et al. (2002). 
At higher temperatures, the spectroscopic Tegf's tend to fall below the structural Tegf's, with the 
reverse situation at lower temperatures and a crossover point at ~1900-2000K. Our derived Teff 
for e Ind B falls nicely on the cool end of this trend. We do not speculate here on the reasons for 
the systematic differences between structural and spectroscopic effective temperatures. 

The treatment of dust in the photosphere is a crucial ingredient in the quantitative spectral 
modelling of the cool L and T dwarfs. Here we have used the unified dust models as discussed 
in detail by Tsuji (2002); however, we also investigated other dust treatments to see what effects 
these would have on the derived physical parameters (primarily Tcgf). Two other sets of model 
atmospheres were generated: 1) one in which dust remained in all layers where the thermochemical 
conditions allowed dust condensation (case B), and 2) the other being the case where all dust 
sank out of the photospheric layers (case C). These very same effects were also investigated and 
discussed in detail by Basri et al. (2000) for a sample of late-M and L dwarfs. In our case B, 
where dust exists over a wide range of depths in the photosphere, the model absorption lines (CO, 
H2O, and CH4) are all much weaker than the observed absorption lines (the addition of significant 
dust opacity over a large region of the photosphere weakens the gas phase absorption lines), and 
no realistic fit to the observed spectra is possible for any reasonable Teff. For the atmospheres 
where all dust sinks from the photosphere (case C), the model molecular absoprtion lines have 
more realistic strengths, but the temperature sensitive CH4 appears at even higher effective 
temperatures; we derive Tcfj=1800K for e Indi Ba. The disagreement with the structural models 
is even worse here. This exercise suggests that some dust in the photosphere provides a better 
physical picture in modelling the high-resolution IR-spectra of e Ind Ba. Our initial comparison 
here points to the need for more high-resolution spectral analyses across the temperature range of 
the L and T dwarfs. 
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4. CONCLUSIONS 

We have used high-resolution infrared spectra of the nearest brown dwarf, e Ind B, to derive 
its physical parameters using comparisons to synthetic spectra calculated from model atmospheres. 
The spectroscopic Teff=1500K, with log g=5.2, and an estimated mass of M=32Mjupiter ■ The 
projected rotational velocity is vsin(i)=28 km s^^, indicating a maximum rotational period of ~3.0 
hours. The comparison between the Tcfj scales derived from spectroscopic plus model atmosphere 
analyses against those derived from Mboi-radius relations reveals a significant systematic difference 
that is unexplained. 

The staff of the Gemini South Observatory are to be thanked for technical support. We 
also thank J. Tennyson for his comments on the H2O assignments, and the referee (G. Basri) 
for helpful suggestions in improving this paper. This research had made use of the SIMBAD 
database, operated at CDS, Strasbourg, France. The work reported here is supported in part by 
the National Science Foundation through AST99-87374 (V.V.S.) and NASA through NAG5-9213 
(V.V.S.). 
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Fig. 1. — The A2.313/iin (top panel) and A1.556/xm (bottom panel) spectra of e Indi B, with the 
plotted wavelengths being those in air. The 2.31//m region is dominated by strong ^^C^^O(2-0) 
lines that are rotationally broadened. The CH4 absorption at A2.318/xm is cither absent or very 
weak at this high spectral resolution. The 1.55//m region exhibits blended features from H2O. 

Fig. 2. — A comparison of the observed spectra of e Ind Ba with synthetic spectra that span a 
range of effective temperatures and gravities. In the top panel the A2.314/im region is shown, with 
the strong ^^C^^O lines, as indicated in Figure 1. In order to fit the observed line shapes the model 
spectra must be broadened by a rotational profile with Vrot=28 km s~^. This spectral region is 
illustrated because of the large temperature sensitivity of the CH4 absorption (as well as the CO). 
A good fit is obtained for Teff=1600K. These particular lines are not very sensitive to gravity over 
the expected values, so a single-gravity set of models is shown (with log g= 5.5). The bottom 
panel shows the observed and model spectra comparisons for a single Tcfj= 1400K but a range of 
gravities. Most of the blended spectral features visible in this wavelength region are from II2O, 
and the relative depth of the absorption is sensitive to gravity, as illustrated. Surface gravities in 
the range of log g= 5.0 to 5.5 are the overall best fits. Note that higher effective temperatures 
will produce extremely weak II2O absorption, while lower Tcfj's result in increasingly dominant 
CH4 absorption, which is not observed. The vertically shifted model spectrum is for Tcff=1300K 
and illustrates the increasingly different absorption as the temperature decreases (note the strong 
feature at 1.5564//m from CH4. 

Fig. 3. — A comparison of effective temperatures derived from model atmosphere spectrum 
synthesis techniques to temperatures derived from structural models for L and T dwarfs: the 
structural Tgfr's are taken from Dahn et al. (2002), while the spectroscopic temperatures are from 
Basri et al. (2000), Leggett et al. (2001), and Schweitzer et al. (2002). There is a well-defined 
systematic trend in the differences between spectroscopic and structural effective temperatures, 
with our result for e Ind Ba falling at the cool end of this trend. 
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